Mon. Not. R. Astron. Soc. 000, 1-?? (2005) Printed 2 February 2008 (MN WT$t style file v2.2) 



Kinematics, ages and metallicities for F and G type stars 
in the solar neighbourhood 

Y. Karatas,, 1 * S. Bilir 1 and W.J. Schuster 2 

I JL " 1 Istanbul University Science Faculty, Department of Astronomy and Space Sciences, 34119 University- Istanbul, Turkey 
, 2 Observatorio Astronomico Nacional, UN AM, Apartado Postal 877, CP. 22800, Ensenada, B.C., Mexico 

o 

1 Accepted 2005 month day. Received year month day; 

ABSTRACT 

A new metallicity distribution and an age-metallicity relation are presented for 
437 nearby F and G turn-off and sub-giant stars selected from radial velocity data of 
Nidever ct al. Photometric metallicities are derived from uvby — H(3 photometry, and 
the stellar ages from the isochrones of Bergbusch & VandenBerg as transformed to 
uvby photometry using the methods of Clem et al. 

The X (stellar-population) criterion of Schuster et al., which combines both kine- 
' matic and metallicity information, provides 22 thick-disk stars, <7w — 32 ± 5 km s _1 , 

I < V rot >= 154 ± 6 km s" 1 and < [M/H] >= -0.55 ± 0.03 dex for these thick-disk 

stars, which is in agreement with values from previous studies of the thick disk, a- 
element abundances which are available for some of these thick-disk stars show the 
typical a-element signatures of the thick disk, supporting the classification procedure 
based on the X criteria. 

Both the scatter in metallicity at a given age and the presence of old, metal-rich 
stars in the age-metallicity relation make it difficult to decide whether or not an age- 
metallicity relation exists for the older thin-disk stars. For ages greater than 3 Gyr, 
our results agree with the other recent studies that there is almost no correlation 
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between age and metallicity, A([M/Fe])/A(age)= -0.01 ± 0.005 dex Gyr -1 . For the 
22 thick-disk stars there is a range in ages of 7-8 Gyr, but again almost no correlation 
between age and metallicity. 

■ For the subset of main-sequence stars with extra-solar planets the age-metallicity 

relation is very similar to that of the total sample, very flat, the main difference being 
that these stars are mostly metal-rich, [M/H]> — 0.2 dex. However, two of these 

stars have [M/H] ~ —0.6 dex and have been classified as thick-disk stars. As for the 
total sample, the range in ages for these stars with extra-solar planetary systems is 
considerable with a nearly uniform distribution over 3 < age < 13 Gyr. 
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1 INTRODUCTION 

Metallicity, kinematic, and age data of F and G type stars 
are of fundamental importance to study the metallicity dis- 
tribution and the age-metallicity relation of the Galactic 
disk. These F and G stars are unique objects for studies of 
the history of the Galactic disk in the sense that photomet- 
ric systems, such as the Stromgren-Crawford uvby — H/3, 
can be used to derive many of their physical parameters, 
such as the interstellar reddening, metallicity, and age. Also, 
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the combined information of stellar photometric abundances 
and kinematics can be used to separate stellar populations 
and to summarize the main properties of a given population. 

With this aim, based on the accurate radial veloc- 
ity observations of Nidever et al. (2002), a data set has 
been created which includes parallaxes, proper motions, and 
uvby — H/3 photometry of late-F, G and K stars. Accurate 
parallaxes and proper motions have become available for 
large numbers of these stars from the Hipparcos and Tycho 
catalogues (ESA 1997) and Tycho-2 (Hog et al. 2000) cat- 
alogue. With this data set, which is large compared to the 
sample of Edvardsson et al. (1993), but small compared to 
the samples of Feltzing, Holmberg & Hurley (2001) and 
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Nordstrom et al. (2004), the metallicity distribution and 
age-metallicity relation in the solar vicinity are determined. 
Much effort has been devoted to the determination of reli- 
able individual isochrone ages for our sample of stars. For 
this task, the isochrones of Bergbusch & VandenBerg (2001) 
and VandenBerg (2003) as transformed to the uvby system 
by the empirically constrained colour-temperature relations 
of Clem et al. (2004), are adopted. 

The sample of Nidever et al. (2002) contains 889 stars 
from the Doppler planet search project at the Keck and 
Lick Observatories (Butler et al. 2000, Vogt et al. 2000). 
This sample contains mostly main-sequence and sub-giant 
stars from F7 to M5, within 50 pc. Chromospherically ac- 
tive stars have been removed from this list as well as stars 
with known stellar companions within 2 " (including known 
spectroscopic binaries). This list provides several advantages 
for our study of the age-metallicity relation: (a) it includes 
many cooler dwarf stars and so avoids one of the selection 
affects discussed by Feltzing et al. (2001); when these cool 
dwarfs are left out of a sample, the old, metal-rich stars are 
preferentially excluded depleting the upper right-hand cor- 
ner of an age-metallicity diagram, (b) the radial velocities 
are especially accurate (systematic errors < ± 0.3 km s _1 ) 
and precise (random errors < ± 0.1 km s _1 ) (Nidever et 
al. 2002) providing us with especially good kinematic data 
for these stars allowing us to separate more cleanly the stel- 
lar populations within the sample, such as the separation 
between the thin- and thick-disk stars, (c) As suggested by 
Gimenez (2000) the uvby photometry can be used to refine 
the survey lists of planetary- search projects and to study 
the metallicities and ages of stars found to have planets. 

This paper is organized as follows: Section 2 describes 
the data; in Section 3 the calculation procedure for the space 
velocities is given; the derivation of astrophysical parame- 
ters: reddening, metallicity, absolute magnitudes, and ages 
is given in Section 4, plus age-metallicity diagrams and a 
discussion of these; and finally the conclusions in Section 5. 



2 DATA 

2.1 Photometry 

For the 889 stars in the catalogue of Nidever et al. (2002), 
uvby photometry of 437 stars and H/3 photometry for 214 
have been collected from the web site of the General Cata- 
logue of Photometric Data 1 of Hauck &: Mermilliod (1998). 
This uvby — Hf3 photometry (b — y, mi, ci, H/3) of these 
437/214 stars are given in Cols. 2-5 of Table 1, respec- 
tively. In Table 1 this uvby photometry of the 437 F and 
G main- sequence, turn-off and sub-giant stars covers the 
ranges 0.087 < mi < 0.497, 0.116 < ci < 0.492, 0.313 < 
(b-y) < 0.580. The 214 stars with H/3 values fall in the 
range, 2.530 < Hf3 < 2.680. Visual magnitudes for our 
sample stars, all brighter than V = 10 mag, and (B — V) 
colours, given in Cols. 7-8 of Table 1, respectively, were 
taken from the Hipparcos and Tycho Reference Catalogues 
(ESA 1997). 



1 http:// obswww.unige.ch/gcpd/gcpd.html 



2.2 Radial Velocities, Parallaxes and Proper 
Motions 

The barycentric radial velocities (V ra d), with a typical accu- 
racy of < 0.3 km s~ , were taken from Nidever et al. (2002). 
Parallaxes (it), the proper motion components (fi a cos5, 
Us) and their associated uncertainties were taken mainly 
from the Hipparcos and Tycho catalogues (ESA 1997) and 
Tycho-2 catalogue (Hog et al. 2000). Parallaxes and their 
uncertainties for three stars (HD 185295, HD 192020 and 
HD 1854) not found in the Hipparcos and Tycho catalogues 
were taken from the catalogue of Kharchenko (2001). The 
radial velocities (V ra d), parallaxes (it), proper motion com- 
ponents (fi a cosS, fis) and their associated uncertainties are 
listed in Table 2. Most of the Hipparcos parallaxes in our 
sample have relative errors that are much less than 20 per 
cent. In Tables 1-3 those stars with extra-solar planetary 
systems have a "P" after their HD number in Col. 1; these 
stars are discussed further in Section 4.4 and especially in 
Fig. 9c. 

2.3 Comparison of our sample with other samples 

Our sample was selected based on the existence of accurate 
radial velocities for 889 FGKM main-sequence and sub-giant 
stars of Nidever et al. (2002). The sample used in this work is 
reduced to only 437 stars, for which uvby— H/3 photometry is 
available, to be able to determine a photometric metallicity 
for each star. Our sample is free from binary contamination 
according to Nidever et al. (2002). The accurate radial ve- 
locities of Nidever et al. (2002) are suitable for the detection 
of extra-solar planets, and so the high accuracy of their ra- 
dial velocities offers us an opportunity to derive very good 
kinematic parameters and thus the age-metallicity relations 
for sub-classes of F and G stars. 

Compared with the samples of Edvardsson et al. (1993) 
and Chen et al. (2003), the size of our sample is large, 
whereas compared to the samples and metallicity distribu- 
tions of Feltzing et al. (2001) and Nordstrom et al. (2004), 
the size of our sample is small but covering a wide metallic- 
ity range, —2.0 <[M/H]< +0.5 dex, with a very high per- 
centage of the stars having —1.0 <[M/H]< +0.5 dex. Our 
sample is not affected by one of the main biases discussed by 
Feltzing et al. (2001) in the sense that old, metal-rich stars 
are amply represented (see Fig. 9 below); the presence of 
these old, metal-rich stars produces an age-metallicity rela- 
tion very similar to other modern age-metallicity relations, 
such as those mentioned above. 

One advantage of having accurate radial velocities for 
our sample together with the ±0.12 — 0.14 dex uncertainty in 
[M/H] (Schuster & Nissen 1989, and see Section 4.3 below) 
is that it allows us to determine well the X cut-offs of the 
stellar populations from Kot and [M/H], where X is the 
stellar-population criterion defined in Schuster, Parrao & 
Contreras-Martmez (1993), a linear combination of V ro t and 
[M/H]. 



3 GALACTIC SPACE VELOCITIES 

Galactic space velocity components (U, V , W) were com- 
puted by applying the algorithms and the transforma- 



Kinematics, ages and metallicities 3 



Table 1. Stellar parameters for our sample stars. The columns arc described in the text. The full version of this table is available online 
at http://www.blackwcllpublishing.com/products/journals/MNR/MNR?/ mnr7sm.htm 



Name 


b - y 


m± 


ci 




E(b-y) 


V 


(B-V) 


M v 


[M/H] 


HD 377 


0.391 


0.204 


0.322 




0.004 


7.59 


0.626 


4.55 


-0.02 


HD 400 


0.332 


0.147 


0.386 


2.615 


-0.014 


6.21 


0.524 


3.61 


-0.35 


HD 1388 


0.382 


0.186 


0.346 




0.001 


6.51 


0.599 


4.41 


-0.08 



Table 2. Hipparcos astrometric, radial velocity and kinematic data of our sample stars. The columns are described in the text. The full 
version of this table is available online at http://www.blackwcllpublishing.com/products/journals/MNR/MNR?/ mnr7sm.htm 



Star Name 


[i a cosS 


Ma 


7T 


v Tad 


U 


V 


W 




(mas yr _1 ) 


(mas yr -1 ) 


(mas) 


(km s" 1 ) 


(km s- 1 ) 


(km b- 1 ) 


(km B- 1 ) 


HD 377 


87.26±0.93 


-0.90±0.53 


25.15±0.97 


1.18 


-14.41 


-7.07 


-3.79 


HD 400 


-114.72±0.53 


-124.60±0.43 


30.26±0.69 


-15.14 


27.64 


-10.36 


-7.85 


HD 1388 


396.97±1.09 


-0.60±0.53 


38.24±0.85 


28.50 


-42.81 


-16.06 


-33.82 



tion matrices of Johnson & Soderblom (1987) to the ba- 
sic data: celestial coordinates (a, 5), proper motion compo- 
nents (nuCosS, us), the radial velocity (Vrad) and the par- 
allax (tt) of each star in Table 2, where the epoch of J2000 
was adopted as described in the International Celestial Ref- 
erence System (ICRS) of the Hipparcos and Tycho Cata- 
logues. The transformation matrices use the notation of a 
right-handed system. Therefore, U, V, W are the compo- 
nents of a velocity vector of a star with respect to the Sun, 
where U is positive toward the Galactic center (I = 0°,b = 
0°); V is positive in the direction of the Galactic rotation 
(/ = 90°, b = 0°); and W is positive toward the north Galac- 
tic pole (b = 90°). U, V, W velocities are heliocentric. The 
velocity components U',V', W have been corrected for the 
local solar motion (U,V,W)® = (+10.0, +14.9, +7.7) km 
s _1 with respect to the local standard of rest (LSR, Allen & 
Santillan 1991) 

The distribution of these Galactic heliocentric velocities 
on the (U, V) plane is displayed in Fig. 1. This Bottlinger 
diagram shows a strong concentration around \U\ < 75 km 
s _1 and —75 < V < 50 km s _1 , where disk-like components 
dominate. However, Fig. 1 also shows more halo-like com- 
ponents having \U\ > 100 km s" 1 and V < — 100 km s _1 . 
Fig. 1 shows that there exist no extreme retrograde veloci- 
ties in our sample; only three stars are retrograde, and these 
all have V > —300 km s _1 (Compare to Fig. 3 of Schuster 
et al. 1993). 
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4 DERIVATION OF ASTROPHYSICAL 
PARAMETERS 

4.1 Reddening 

The interstellar extinction for our sample stars has been 
measured due to its effect on the derivation of the metal- 
licities and absolute magnitudes. E(b — y) of 214 stars with 
2.530 < Hf3 < 2.680 were estimated from an intrinsic colour 
calibration of (b — y) vs H/3 (Schuster & Nissen 1989). A 
few iterations were made to obtain E(b — y), with a consis- 
tency of ±0.0001 for (6 - y) . For 223 stars which lack H/3, 
E{B — V)s was obtained from the Schlegel, Finkbeiner & 
Davis (1998) reddening maps. We adopt a V band absorp- 
tion Ay — 3.1E(B — V) and assume that the dust layer has 



Figure 1. The Bottlinger diagram for our sample. The velocities 
arc heliocentric. 



a scale height H — 125 pc. Thus, the reddening E(B — V)s 
for a given star is reduced compared to the total reddening 
E(B - V) T by a factor (1 - e -^in(b)/H^ where b &nd d afe 

the Galactic latitude and distance, respectively. Note, how- 
ever, that Arce & Goodman (1999) caution that Schlegel et 
al. (1998) reddening maps overestimate the reddening val- 
ues when the colour excess E(B — V)s is more than 0.15 mag. 
Hence, we have adopted (Schuster et al. 2004) a slight revi- 
sion of the Schlegel et al. (1998) reddening estimates, via an 
equation, E{B-V)a = 0.10 + 0.65(S(B - V)s - 0.10) when 
E(B-V)s > 0.10, otherwise E(B-V) A = E(B-V)s, where 
E(B — V)a indicates the adopted reddening estimate. Figure 
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Figure 3. Differences of the mctallicities derived from the for- 
mulae of Martell & Laughlin (2002) minus those from Schuster 
& Nissen (1989) versus the values from Schuster & Nissen (1989) 
([M / H]g]<i) , for the stars with (b — y)o > 0.47 in our sample. 
Note that the systematic errors become obvious only for [M / -ff]sN 
greater than about solar, maximize at [M/if]gN ~ +0.07 dex, and 
then become small again at [M/H]sn ~ +0.15 dex. 



magnitudes of our sample are presented in Col. 9 of Table 
1. 



Figure 2. E(b — y) histogram of all stars in our sample. 



3 of Schuster et al. (2004) shows that there is no significant 
systematic difference between the reddenings obtained from 
these two methods, the intrinsic colour calibration and the 
reddening maps. The final adopted E(B — V)s values of 
223 stars are transformed into E(b — y) via the equation 
E(B -V) = 1.35£(6 - y) given by Crawford (1975). 

Fig. 2 shows the E(b — y) distribution of the 437 stars. 
As expected, the stars with a median distance of 36 pc show 
a strong peak at about E(b — y) = 0.000, with somewhat 
asymmetric wings, the negative tail being slightly bigger 
than the positive, displaced only 0.005 mag from symmetry. 
Reddenings less than 0.015 are mostly due to the observa- 
tional errors of the photometry (Nissen 1994). For our final 
sample, reddening corrections have been applied to only six 
stars, those with E(b — y) > 0.015; for all the remaining 
stars no reddening correction has been applied. For these 
six stars, the indices ci and mi are dereddened by the re- 
lations, c = ci — 0.2E(b — y) and m = mi ± 0.3E(b — y), 
given by Crawford (1975). The final reddenings, E(b — y), 
are given in Col. 6 of Table 1. 

4.2 Absolute magnitudes 

The absolute visual magnitudes, My, are needed to deter- 
mine the ages of stars in our sample from the Mv, (b — y) 
plane. Although relative errors (a^/n) of Hipparcos paral- 
laxes of the stars in Table 2 are less than 0.20, the Lutz- 
Kclkcr correction (Lutz & Kelker 1973) has been applied to 
the absolute magnitudes of our sample using the following 
relation constructed from their data: 

SM V = a(a n /n) 3 + b(a n /n) 2 + c(ov/tt) + d (1) 

where a = -63.0303, b = -1.0736, c = -0.3286, d = 0.0008. 

The relation M v = V ~5log(l/n)+5+8M v allows us to 
determine absolute magnitudes when Hipparcos parallaxes 
of the stars in our sample are available. Here, V magnitude 
is corrected for reddening using V = V — Ay- The absolute 



4.3 The metallicity distribution and the 
separation of population types 

In order to estimate the metallicities of our sample stars 
from the uvby — HP photometry presented in Table 1, the 
metallicity ([Fe/H]) calibrations of Schuster & Nissen (1989) 
have been used. We assume here that [M/H] = [Fe/H]. 
However, Twarog, Anthony- Twarog & Tanner (2002) re- 
port that the metallicity calibration of G stars given by 
Schuster & Nissen (1989) has a systematic error that is 
strongly colour dependent, reducing the estimated [M/H] 
for stars above [M/H] ~ —0.2 dex for the very reddest 
stars, (6 - y) > 0.47. Martell & Laughlin (2002) give an 
improved uvby — H/3 metallicity relation for F, G and early 
K stars. In order to see this systematic effect, in Fig. 3 the 
differences between the metallicities derived from Martell & 
Laughlin (2002) minus those from Schuster & Nissen (1989) 
have been plotted versus those derived from the Schuster 
& Nissen (1989) formulae. In Fig. 3 it is worth noting that 
only the stars with (b - y) > 0.47 and [M/H] > - 0.03 dex, 
show a significant difference, increasing to a maximum at 
[Fe/H] w +0.07 dex and then decreasing again. So, for the 
final metallicities of this paper, the [Fe/H] calibration given 
by Martell & Laughlin (2002) has been applied to only the 
stars with (b - y) > 0.47 and [M/H] > -0.2 dex; for the 
remaining stars, the metallicities estimated via the formulae 
of Schuster & Nissen (1989) are adopted. The metallicities 
of our stars are given in the Col. 10 of Table 1. 

In Fig. 4, for 71 stars, our photometrically derived 
[M/H] values are plotted against the spectroscopic values of 
the Cayrel de Strobel, Soubrian & Ralite (2001) catalogue. 
A least square fit to the points of Fig. 4 yields a regression 
curve of [M/H] spec = 0.99(±0.04)[M/ii"] phot - 0.01 (±0.02) 
for the metallicity range -2.0 <[M/H]< +0.5 dex. The 
agreement is excellent with a scatter of a = ±0.12, exclud- 
ing the large deviation of HD 202573, which is a giant star 
(Henry et al. 2000) and which is outside the range of the 
Schuster & Nissen (1989) [Fe/H] calibrations.. It is clear 
from Fig. 4 that accurate metal abundances can be derived 
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Figure 4. uvby — H(3 photometric metallicities of 71 stars versus 
the spectroscopic ones taken from the Cayrel de Strobel ct al. 
(2001) catalogue. The one-to-one relation is drawn as a solid line. 

from uvby — H f3 photometry, using the combined calibrations 
as described above. 

For the empirical [Fe/H] calibrations of Schuster & Nis- 
sen (1989) the estimated standard deviations of a single 
photometric determination of [Fe/H] were ±0.14 at [Fe/H] 
« -0.5 dex, and ±0.21 at [Fe/H] « -1.5 dex. Compar- 
isons by Feltzing et al. (2001) of photometric abundances 
from the Schuster & Nissen (1989) calibration equations 
with abundances from two recent spectroscopic studies Ed- 
vardsson et al. (1993) and Chen et al. (2000) have shown 
that these error estimates are overly conservative; they find 
a scatter of only ±0.10—0.11 for their more metal- rich group, 
in very good agreement with the value of ±0.12 found above 
from a somewhat less homogeneous spectroscopic data set. 

In Fig. 5 the distribution of metallicity, [M/H], for our 
sample stars is presented, and also shown is a fit to the his- 
togram using three Gaussians. The sample is mix of all the 
stellar populations that are represented in the solar neigh- 
bourhood, the thin and thick disks, and the halo. As is ex- 
pected, in the solar neighbourhood, the majority of the stars 
belong to the disk with the thin disk dominating, a signifi- 
cant but smaller contribution from the thick disk, and only 
a few halo stars (~ 8; see Fig. 6 below). In Fig. 5 two main 
components can be seen, one with —0.9 < [M/H] < ±0.5 dex 
(the disk contribution: thin plus thick) and the other with 
[M/H] < — 0.9 dex (mostly the halo stars). A rough Gaus- 
sian fit has been made to this metallicity histogram using 
the mathematical package 'Origin' 2 . One disk-like compo- 
nent has < [M/H] >~ —0.04 dex, a dispersion of 0.37, and 
is nearly symmetric. Our sample contains an obvious contri- 
bution from the thick disk with a Gaussian having a mean 
abundance of < [M/H] >~ —0.45 dex and a dispersion of 
0.41. There are eight halo stars with < [M/H] > < - 0.9 

2 http:// www.originlab.com 
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Figure 5. Distribution of the metallicities of our stars. Curves 
are the fitted Gaussian distributions, which give the mean metal 
abundances of < [M/H] >= -0.04, -0.45, and -1.30 dex for 
the thin and thick disks, and halo (thin curves), respectively, and 
their sum (heavy curve). 

dex (see also Fig. 6). Nevertheless, thin and thick-disk stars 
are not clearly distinguishable from Fig. 5. There exists no 
clean, straightforward procedure for separating the stars of 
these different disk populations, which overlap considerably 
in the metallicity distribution. Kinematics are usually in- 
voked to improve the thin/thick disk separation, and more 
recently [a/Fe] abundances (see for example, Bensby, Feltz- 
ing & Lundstrom 2003). 

In Fig. 6 the Kot, [M/H] diagram for our sample stars is 
plotted. Two components are clearly noted: the halo compo- 
nent, Kot < 100 km s _1 and [M/H] < - 0.9 dex; the other 
disk-like, centered at Kot ~ 215 km s _1 and [M/H] ~ 0.0 
dex. The rotation velocity of the LSR about the Galactic 
center is taken here to be 220 km s _1 so that Kot = V'+ 220 
km s _1 is the rest-frame rotation velocity of a given star 
(Kerr & LyndenBell 1986). More modern values for the cir- 
cular speed of the Milky Way at the solar circle are some- 
what larger, such as Oo = 234±13 km s _1 given by Fukugita 
& Peebles (2004). Here, V is the velocity with respect to 
the LSR. It can be seen from Fig. 6 that the majority of 
the stars have Kot ~ 200 km s _1 and [M/H] > -0.5 dex. 
These stars are thin-disk stars. Even though there is scarcity 
of stars with [M/H] < -0.5 dex and Kot < 150 km s -1 , 
Kot increases approximately linearly with [M/H] between 
-1.0 < [M/H] < -0.5 dex. For the [M/H] < -1.0 dex, it 
appears that there is no correlation between [M/H] and Kot, 
but there are too few stars for any firm conclusion. 

This Kot, [M/H] diagram can be used to separate out 
the different stellar populations, as discussed in Nissen & 
Schuster (1991) and in Schuster et al. (1993). In the former 
reference, a diagonal cut connecting ([M/H], Kot) = (—0.3, 
km s _1 ) and ( — 1.5, 175 km s _1 ) was used to separate halo 
stars from 'high-velocity disk' stars. In the latter reference, 
the parameter "X" was defined, which is a linear combina- 
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Table 3. Kinematic and abundance data of thick-disk stars. The parameter X defined by Schuster ct al. (1993) in Col. 3 denotes the 
linear combination of V ro t and [M/H] prlo t- [Fe/ H] Bpcc , a-elements and [a/Fe] are given in Cols. 6-10. [et/Fe] is defined as (l/4)([Mg/Fc] 
+ [Si/Fc] + [Ca/Fc] + [Ti/Fe]). The references are as follows: 1, Bcnsby et al. (2005); 2, Pompeia, Barbuy & Grenon (2003); 3, Cenarro 
et al. (2001); 4, Cayrel de Strobel et al. (2001); 5, Fulbright (2000); 6, Mishcnina et al. (2004); 7, Fricl (1987); 8, Santos, Israelian & 
Mayor (2004) 



Name Vrot X [M/H] phot [Fe/H] spcc [Mg/Fe] [Si/Fc] [Ca/Fe] [Ti/Fc] [a/Fe] Refs. 

(km s- 1 ) 



HD 3795 


145 


-18 


-0.39 


-0.59 


0.31 


0.22 


0.24 


0.29 


0.27 


1 


HD 16623 


143 


-16 


-0.50 


-0.60 


0.25 


0.10 


-0.02 


0.00 


0.08 


2 


HD 19034 


166 


-20 


-0.41 
















HD 22879 


149 


-12 


-0.74 


-0.84 


0.40 


0.27 


0.27 


0.31 


0.31 


1 


HD 29587 


185 


-20 


-0.59 


-0.61 


0.31 


0.17 


0.17 


0.19 


0.21 


1 


HD 30649 


155 


-16 


-0.56 


-0.50 






0.08 






3 


HD 32387 
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tion of Vrot and [M/H], and is used to make other diagonal 
cuts which isolate more cleanly the thick-disk stars. Our 
most recent work (Schuster et al. 2005) indicates that the 
range —21 < X < —6 gives a fairly clean thick-disk sample, 
with only small contamination by the halo and old thin disk. 
In Schuster et al. (1993) the cut -21 < X < -18 was used 
to define an even cleaner thick-disk sample, but here too few 
stars are found in this reduced X interval. 

The X criterion provides a cleaner separation of the 
thin- and thick-disk components even though they overlap 
considerably in both Kot and [Fe/H]. Using the definition 
of X given by Schuster et al. (1993), combining Kot and 
[Fe/H], the X values of the stars have been calculated. The 
X distribution of the stars is given in Fig. 7. 

Comparing our X distribution with the one plotted by 
Schuster et al. (1993), constructed from a larger data base, 
X cut-offs for the three stellar components are fixed for 
Fig. 7. The number of stars in the thick-disk sample, de- 
fined by the range, —21 < X < —6, is 22. These limits for 
the thick disk have been derived from the work of Schuster 
et al. (1993) and from a project that is in process (Schuster 
et al. 2005). The thin disk dominates for X < —33, and the 
halo with a small fraction of the total extends for X > —5. 
These fixed cut-offs of three stellar populations are indicated 
as dashed vertical lines on Fig. 7. 

The vertical velocity dispersion, aw = 32 ± 5 km s _1 
(error calculated according to Appendix E of Taylor 1997), 
is estimated for the range —21 < X < —6, where the thick 
disk dominates. This value is less than the ctw = 46 km 
s _1 of Schuster et al. (1993), but is in agreement with 
the range of 30-37 km s _1 given by Norris (1987), Croswell 



et al. (1991), and Carney, Latham & Laird (1989). Other 
values, < Vrot >= 154 ±6 km s _1 and < [M/H] >= 
—0.55 ± 0.03 dex (mean errors) for this thick-disk sample 
are in concordance with the mean values given in the liter- 
ature, such as in Schuster et al. (1993). 

Our thick-disk sample defined by the range, —21 < X < 
—6, is given in Table 3 together with the kinematic and 
abundance data. The a-element and spectroscopic [Fc/H] 
abundances of some stars have been found in the literature. 
It can be seen from Table 3 that spectroscopic [Fe/H] abun- 
dances are in good agreement with the photometric ones. 
The mean < [a/Fe] > values of seven thick-disk stars, ex- 
cept HD 16623 are constant at a level of 0.27 dex, which are 
overabundant relative to thin disk stars, over the spectro- 
scopic metallicity range from -0.45 to -0.84 dex (See Fig. 1.4 
in Nissen 2004, Fig. 4 in Nissen & Schuster 1997, Fig. 13 in 
Bensby et al. 2003, and Fig. 8 in Bensby et al. 2005). 

[Mg/Fe] abundances of stars with thick-disk kinemat- 
ics in Table 3, which fall in the range +0.19 < [Mg/Fe] < 
+0.43 dex are similar to those of the thick disk found in 
the literature. The trend of Q-capture elements like Mg, 
Si, Ca, and Ti with photometric [M/H] (or spectroscopic 
[Fe/H]) abundances of thick-disk stars in Table 3 agrees well 
with the findings of Nissen & Schuster (1997) and Bensby 
et al. (2005). Also, the a-element abundances of thick-disk 
stars, which are available in Table 3, support the classifica- 
tion procedure based on X criteria defined by Schuster et 
al. (1993). 

The only expection might be HD 16623; its a/Fe is 
low by a factor of 1.5 relative to the seven other thick-disk 
stars in Table 3. The kinematics and photometric (and spec- 
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and thick disks, and the halo. 



troscopic) [M/H] abundance of HD 16623 indicate that this 
star is a a member of the thick disk. The mean [a/Fe] of 
HD 16623 is low due to its [Ca/Fe] and [Ti/Fe] abundances, 
which are more similar to that of the thin disk, whereas 
its [Mg/Fe] abundance is higher and more similar to that 
of the thick disk. This star is perhaps transitional between 
the thick and thin disks, with a more mixed chemical and 
dynamical history. 

Another star which stands out in Table 3, and also in 
Fig. 6, is HD 104556. Its X, [M/H] phot , and W (= +15 km 
s _1 , from Table 2) values would all indicate a (thick) disk 
star while its V TO t value (+74 km s _1 ) is more indicative of 
the halo. For example, from Table 5 of Schuster et al. (1993) 
the Kot value of HD 104556 is 2.1a removed from the mean 
value of the thick disk, but only 0.5ct from that of the halo. 
Other studies have also given compositions for HD 104556 
which are more thick-disk-like than halo-like, such as Friel 
(1987), who gives [Fe/H] = -0.64 dex, and Eggen (1997), 
who gives [Fe/H] = —0.55 dex. These stars point out the 
risks in using only a single parameter, such as X, [Fe/H], or 
[a/Fe], to derive the population type of an individual star, 
and also show that some stars do not fit cleanly into any 
one population type, according to all criteria. 

4.4 Ages and an age-metallicity relation 

Fig. 8 shows the position of all the thin-disk stars in our sam- 
ple, with a mean metal abundance of < [M/H] >= —0.04 
dex for the metallicity range -0.43 < [M/H] < +0.35 dex, 
as compared with the isochrones of VandenBerg (2003) as 
transformed to the uvby photometry by Clem et al. (2004). 
It is seen that the majority of the stars are distributed in- 
side the isochrones. Possible systematic shifts between the 
isochrones and the positions of the redder main-sequence 
stars have been examined, and no need for systematic cor- 
rections was found. It is also conspicuous that a few stars 
with the values of (6 - y) D >0.44 and M v >4.5 lie off the 
isochrone grid, i.e. ages > 16 Gyr. These may be unidentified 
binary stars. 

Isochrone ages are determined by placing the stars in 
observational HR diagrams (such as Fig. 8) using the ob- 
served My, (b — y) values for each star, and then reading 
off the age of the star by interpolating between the theoret- 
ically computed isochrones. The isochrones of VandenBerg 
(2003) as transformed to the uvby system by the colour- 
temperature relations of Clem et al. (2004) have been cho- 
sen for the age determination of the thin-disk stars. These 
uvby isochrones cover the range —0.60 < [M/H] < +0.50 
dex and ages from 1 to 20 Gyr in 1 Gyr steps. Whereas, 
for the age determination of the thick-disk and halo stars, 
which are separated according to the X criterion, the set 
of isochrones from Bergbusch & VandenBerg (2001), as 
transformed by the same colour-temperature relations, have 
been used. This latter set of isochrones cover —2.31 < 
[M/H] < -0.30 dex and ages from 6-20 Gyr in 1 Gyr steps. 
The ages of the 22 thick-disk stars have been interpolated 
in the My, (b — y) diagram using a constant [a/Fe] = 
+0.30 dex and interpolated for the metal abundances of 
[M/H] = -0.30,-0.40,-0.52,-0.61,-0.70,-0.83,-1.01 
dex. At sub-solar [M/H] we assumed [a/Fe] = +0.30 dex 
for thick-disk stars with [M/H] < — 0.8 dex, taking into 
consideration Fig. 3 of Wheeler, Sneden & Truran (1989). 
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In order to estimate the errors for the age interpolations, 
the uncertainties (cry, and (Jb-y) in V, and (b — y) have been 
collected from the web site of the General Catalogue of Pho- 
tometric Data of Hauck & Mermilliod (1998), and from the 
Hipparcos and Tycho catalogues (ESA 1997), the uncer- 
tainties (a n /n) in the parallaxes (%) of our sample of stars 
given in Table 2. The derived error equation for the Pogson 
relation given in Section 4.2, om v = cry + 2.17(cr x /7r), al- 
lows us to estimate ctm v from the uncertainties in V and tt of 
individual stars in our sample. The resulting average uncer- 
tainties in My and (b — y) are 0.081 and 0.003 mag, respec- 
tively. For stars lying along the turn-offs of the isochrones, 
these average uncertainties in My and (b — y) result in age 
errors of ~ 1 Gyr when applied to the difference between 
two isochrones in the My, (b — y) diagram. For sub-giants 
above the turn-offs, the errors are somewhat larger, 1.0 — 1.5 
Gyr. Below the turn-offs the age errors can grow indefinitely 
due to the convergence of the isochrones. We have measured 
ages to an error limit of about ±2.5 Gyr. Error bars corre- 
sponding to an uncertainties of 1 Gyr in age and 0.12 — 0.14 
dex in [Fe/H], (Schuster & Nissen 1989 and Section 4.3) are 
indicated in Fig. 9a. 

In Fig. 9a the metallicity is plotted as a function of age 
with different symbols for the definitive thin disk (circles, 
X < —33), probable thin disk (cross signs, —33 < X < —21) 
and probable thick disk (plus signs, —21 < X < —6). It can 
be seen from Fig. 9a that few stars have ages less than 3 
Gyr; this is to be expected since the stars of this sample by 
Nidever et al. (2002) contain mostly late F-, G-, and K-type 
stars. This derived age-metallicity relation agrees qualita- 
tively with that of Edvardsson et al. (1993). 



Figure 9. (a) Age-metallicity relation for our stellar sample. The 
circles show the stars with X < —33, which arc the definite thin- 
disk stars, Cross signs indicate stars with —33 < X < —21, which 
are probable thin-disk stars but with some small probability of 
being thick disk, and the plus signs show the probable thick-disk 
stars, respectively. A typical average error bar, referring to the 
uncertainties of 1 Gyr in age and 0.12 — 0.14 dex in Schuster & 
Nissen (1989)'s metallicity calibration, is indicated, (b) Average 
metallicity < [M/H] > in age bins. The straight line is a least- 
squares fit having a slope of —0.01 ± 0.005 dex Gyr -1 . Error 
bars are the standard deviation of the mean of the values in each 
bin. (c) Age-metallicity relation for 38 main-sequence stars which 
have extra-solar planets according to the catalogues of Schneider 
(2005) and Mayor et al. (2005) 



A slight trend that metallicity decreases with increasing 
age is seen in Fig. 9b; the slope of a least squares fit to the 
< [M/H] > versus age is -0.01 ± 0.005 dex Gyr -1 . Even 
though there is this overall trend of decreasing mean metal- 
licity with increasing age, both a larger scatter in [M/H] at 
a given age and the presence of old, metal-rich stars make it 
difficult to decide whether or not an age-metallicity relation 
really exists for the thin disk. So, it is difficult to draw any 
clear conclusion from Fig. 9a. 

The age-metallicity relations derived from the larger 
samples of Feltzing et al. (2001) and Nordstrom et al. (2004) 
apply mainly for stars less than 3 Gyr, whereas ages of our 
stars are generally greater than this limit. For ages larger 
than 3 Gyr the age-metallicity diagrams of Feltzing et al. 
(2001; their Fig. 10) and of Nordstrom et al. (2004; their 
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Fig. 27) appear very similar to our Fig. 9a, almost no cor- 
relation between age and metallicity. 

In Fig. 9a the 22 thick-disk stars show a range of 7- 
8 Gyr in their ages but no obvious correlation between the 
ages and metallicities, in contrast with the results of Bensby, 
Feltzing & Lundstrom (2004), who find a pronounced age- 
metallicity relation for the thick disk. However, Bensby et 
al. (2004) work with a thick-disk sample ten times larger 
than ours, separated from the thin disk more rigorously, and 
so their conclusions should be preferred. We do agree with 
them that the thick-disk stars show a very significant range 
in ages, > 5 Gyr. 

From Fig. 9a, it can be seen that there is lack of metal- 
poor young stars in the left-bottom corner. However, there 
exists numerous old metal-rich stars in the right-upper cor- 
ner, which leads to the substantial scatter and the flatten- 
ing of the plot. The existence of the old metal-rich stars can 
also seen in the works of Feltzing et al. (2001), Nordstrom et 
al. (2004), Rocha-Pinto et al. (2000), and Chen et al. (2003). 

Fig. 9c shows the same age-metallicity plot for the main 
sequence stars from Nidever et al. (2002) for which extra- 
solar planets have been detected (Schneider 2005; Mayor et 
al. 2005). Again, the relation between the metallicity and 
the age is very flat as in Fig. 9a, the main difference being 
that the large majority of these stars have high metallici- 
ties, [M/H] > — 0.2 dex, as has been pointed out in many 
studies, such as Gonzalez (1997), Laughlin (2000), Santos, 
Israelian & Mayor (2001), Murray & Chaboyer (2002), and 
Reid (2002). However, two of these stars (HD 114729 and 
HD 114762) have lower metallicities, [M/H] ~ -0.6 dex, 
and have been classified thick disk, with rather large ages, 
~ 12 Gyr. Interestingly, this entire sample of stars with 
extra-solar planetary systems has large ages and a large 
range in ages, 3 < age < 13 Gyr. 

In Fig. 10 we plot the ages and metallicities for our total 
sample, for three bins of My. A similar pattern as in Fig. 14 
of Feltzing et al. (2001) appears. As we allow stars with faint 
absolute magnitudes to enter our samples, the upper right- 
hand corner of the age-metallicity plot is progressively filled 
in. In panel (c) of Fig. 10, stars with [M/H] > -0.1 dex 
and ages larger than 10 Gyr become numerous in the up- 
per right-hand corner. From panel (b) of Fig. 10, one might 
conclude a fairly obvious and pronounced age-metallicity re- 
lation, but with the stars of panel (c) this has mostly disap- 
peared. Fig. 10c reveals that old, metal-rich stars are domi- 
nant at the faint absolute magnitudes. 

Chen et al. (2003) investigate the nature of these old 
metal-rich stars in the solar neighbourhood and conclude 
that most seem to have originated from the inner thin disk. 
They note that HD 190360 is an exception in their sample 
and suggest that it is a metal-rich thick-disk star. A 15.8 
Gyr age can be attributed to this star since its position in 
the Mv, (b—y) diagram of Fig. 8 places it slightly to the left 
of the 16 Gyr isochrone. Our photometric metal abundance 
([M/H] = +0.17 dex) and a 15.8 Gyr age show that this star 
is old and metal-rich, in accordance with the suggestion of 
Chen et al. (2003). For this star, our kinematics (V ro t =190 
km s _1 , W' — —56 km s _1 ) and photometric abundance 
([M/H] — +0.17 dex) are almost identical to the ones of 
Chen et al. (2003). Whereas our X — —35 value indicates 
that it is a thin-disk star, the rather large absolute value for 
W' would suggest thick disk, again in agreement with Chen 
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Figure 10. Age-metallicity plots of three sub-samples for the (a) 
stars with My < 3, (b) stars with 3 < My < 4, (c) stars with 
4 < My < 6.3. The symbols are the same as in Fig. 9a. 



et al. (2003). HD 190360 is another star which does not fit 
cleanly into all the criteria for a single stellar population, 
and shows again the risk in using only two variables, [M/H] 
and Kot, to define the population type of a single star. 

The scatter observed in age-metallicity plot of Fig. 9a 
can be interpreted within the scope of models of Galactic 
chemical evolution of Pilyugin & Edmunds (1996) and Rai- 
teri, Villata & Navarro (1996). These models consider the 
chemical enrichment with time, together with the dynam- 
ical evolution. The model of Pilyugin & Edmunds (1996) 
specifically takes into account both the self-enrichment of 
star forming regions with sequential star formation and the 
effect of irregular rates of infall of un-processed matter onto 
the stellar disk. These Galactic chemical-evolution mod- 
els can generate a very conspicuous scatter, such as that 
seen in the observational age-metallicity plots, such as our 
Fig. 9a, Fig. 14a of Edvardsson et al. (1993), Figs. 27-28 
of Nordstrom et al. (2004), and Figs. 10-11 of Feltzing et 
al. (2001). However, for the models of Pilyugin & Edmunds 
(1996) episodic gas infall onto the disk is required to pro- 
duce a scatter in metallicity ([Fe/H]) as large as that ob- 
served; self-enrichment processes alone will not do it. The 
N-body/hydrodynamical simulations of the chemical evolu- 
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tion by Raiteri et al. (1996) can produce a scatter in [Fe/H] 
even larger than that observed (see their Fig. 9) depending 
upon the volume in the Galaxy over which the results have 
been considered. 

Our sample based on radial velocities of Nidever et al. 
(2002) arc not well populated at all ages, but contains late- 
F, G and K stars with apparent magnitudes fainter than 
V — 8 mag, as compared to the samples of Edvardsson et 
al. (1993), Feltzing et al. (2001) and Nordstrom et al. (2004). 
This selection effect causes a bias in the sense that ages less 
than 3 Gyr are lacking, but the sample is without kinematic 
bias. Moreover, the metallicity distribution of our sample is 
representative of the older stellar populations of the solar 
neighbourhood, and shows a small overall trend in the age- 
metallicity plot. The presence of old, metal-rich stars, from 
the cooler program stars, shows that this sample has avoided 
one of the distortions of the age-metallicity relation, such as 
that which affected the results of Twarog (1980), as pointed 
out by Feltzing et al. (2001). Also, the presence of a signif- 
icant number of metal-rich and older stars emphasizes the 
usefulness of the planetary-search list of Nidever et al. (2002) 
since many studies, such as those of Gonzalez (1997), Laugh- 
lin (2000), Santos et al. (2001), Murray & Chaboyer (2002), 
and Reid (2002), have shown that the parent stars of most 
extra-solar planets are significantly metal-rich, as confirmed 
in our Fig. 9c. 



5 CONCLUSIONS 

Our main conclusions are as follows: 

1) Accurate metal abundances can be derived from 
the uvby — Hp photometry, using the combined calibra- 
tions as described above in Section 4.3. In Fig. 4 a com- 
parison between the photometric metallicities and modern 
spectroscopic ones for 71 stars shows excellent agreement: 
[M/H] spcc = 0.99(±0.04)[M/ii"] phot - 0.01(±0.02) for the 
metallicity range —2.0 < [M/H] < +0.5 dex with a scatter 
of ff[M/H] = ±0.12. This error agrees well with estimates 
made by Feltzing et al. (2001), ±0.10-0.11. 

2) The metallicity and V TO %, [M/H] distributions of our 
sample show a mixture of all the stellar populations in the 
solar neighbourhood; the thin disk dominates with a small 
but significant contribution from the thick disk, plus a few 
halo stars. 

3) Using the X criterion of Schuster et al. (1993) has 
allowed us to isolate a thick-disk sample of 22 stars with 
minimal contamination from the other stellar populations. 
A aw = 32 ± 5 km s _1 velocity dispersion is estimated 
for the range —21 < X < —6 where these thick-disk stars 
dominate. This value is within the range of 30-37 km s _1 
given by Norris (1987), Croswell et al. (1991), and Carney 
et al. (1989). These 22 thick-disk stars also have < Kot >= 
154 ± 6 km s" 1 and < [M/H] >= -0.55 ± 0.03 dex. 

4) a-clcment abundances which are available for some 
thick-disk stars in our sample show the typical tv-element 
signatures of the thick disk, supporting the classification 
procedure based on the X criteria, which is a linear com- 
bination of Vrot and [Fe/H] (Schuster et al. 1993). 

5) There is a lack of stars with ages less than 3 Gyr 
in our age-metallicity diagram, due to the selection criteria 
of the sample. For ages greater than 3 Gyr, both a larger 



scatter in [M/H] at a given age and the presence of old, 
metal-rich stars make it difficult to decide whether or not 
an age-metallicity relation really exists, very similar to the 
works of Edvardsson et al. (1993), Feltzing et al. (2001) and 
Nordstrom et al. (2004). Our age-metallicity plot (Fig. 9) 
shows a very shallow slope of A([M / Fe]) / A(age) — —0.01 ± 
0.005 dex Gyr^ 1 . 

6) Our Fig. 9a shows a significant age span for the thick 
disk with a range in ages > 5 Gyr, in agreement with the 
work of Bensby et al. (2004). 

7) The age-metallcity diagram for the stars with extra- 
solar planets shows a relation very similar to that of the 
total sample, very flat with almost no variation in the mean 
metallicity with age; the main difference is that most of these 
stars are metal-rich, [M/H] > —0.2 dex with the exception of 
two (HD 114729 and HD 114762) which have [M/H] ~ -0.6 
dex and have been classified as thick disk. This subsamplc 
also shows a more or less uniform distribution over a wide 
range in ages, 3 < age < 13 Gyr, similar to the total sample. 

8) A number of stars have been pointed out which do 
not fit well into any stellar-population classification scheme, 
such as HD 16623, HD 104556, and HD 190360. These stars 
exemplify the risks in using only a single criterion, such as X, 
[Fe/H], Kot, [a/Fe], or W' , to classify the stellar-population- 
type of any individual star. All such information should be 
used, but even then, some stars cannot be classified unam- 
biguously. 
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